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Abstract 

We use the parabolic wave equation to study the propagation of x-rays in tapered waveguides by numercial simulation 
and optimization. The goal of the study is to elucidate how beam concentration can be best achieved in x-ray 
optical nanostructures. Such optimized waveguides can e.g. be used to investigate single biomolecules. Here, we 
compare tapering geometries, which can be parametrized by linear and third-order (Bezier-type) functions and can 
be fabricated using standard e-beam litography units. These geometries can be described in two and four-dimensional 
parameter spaces, respectively. In both geometries, we observe a rugged structure of the optimization problem's "gain 
landscape" . Thus, the optimization of x-ray nanostructures in general will be a highly nontrivial optimization problem. 

1. Introduction 

X-ray waveguides are promising novel optical devices to generate x-ray quasi-point sources for hard x-ray 
imaging |1|2] , X-ray waveguides (WG) have been designed and fabricated as planar layered systems (1D-WG) 
[3I4I5| or as lithographic channels (2D-WG) 6 7. In both cases (guiding layer or guiding channel) mainly 
waveguides with a core of constant cross section have been considered. Beam propagation in such structures, 
which are translationally invariant along the propagation axis z, is by now well studied and understood 
[819110111] . As optical elements, such waveguides serve mainly two purposes: Firstly, the radiation modes 
are damped out and the cross section of the exit beam is thus reduced significantly with respect to the 
incident beam. Secondly, in the case of single- mode waveguides |12j . the phase and amplitude of the wave is 
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determined exactly, and therfore presents an ideal wave for coherent imaging. In other words, the waveguide 
acts as a spatial and coherence filter, but not to collimate or to concentrate the beam. The flux density 
necessary for nanometer-sized beam with high enough intensities for practical applications must then be 
achieved by focusing elements placed in front of the waveguide, such as Kirkpatrik-Baez (KB) mirrors |13|7j . 
Fresnel zone plates [12] or compound refractive lenses [15116] . 

Here we address the question to which extent tapering of the waveguide can be used to concentrate the 
beam and to filter it at the same time. This would make pre-focusing optics in front of the waveguide obsolete 
for some applications. A first experimental realization was demonstrated by tapered air-gap waveguides 
[17j . but was restricted to the multi- modal regime. Bergemann et al. have then studied linearly tapered 
x-ray waveguides by analytical and numerical calculation [18 , concentrating on minimum beam width at 
the exit. Here we want to generalize the numerical studies to non-linear tapering profiles. Rather than 
the minimum achievable beam width, which has been addressed before [18] , we address the maximum flux 
density enhancement as function of tapering parameters, and the associated question of an optimum interface 
geometry. As in all problems of reflective optics, e.g. bent mirrors, one would assume that the shape function 
is a crucial parameter in the problem. As a first step, we consider both linear and third order polynominal 
interface geometries. Note that the parameters studied correspond well to structures which can also be 
fabricated by e-beam lithography. The propagation in the stuctures is studied by numerical solution of the 
parabolic wave equation (PWE), as used previously for standard waveguides [19 10J. The use of the PWE 
will be briefly explained in the next section, followed by the simulation results. 

2. The simulation method 

The propagation of monochromatic x-rays in materials is described by the Helmholtz equation 

[A + k 2 n 2 ]^ = . (1) 

'J denotes an electric field component in this scalar wave equation, k = ^ the wavenumber, and A the 
wavelength. The index of refraction is given by n = 1 — 5 — i/3, where the imaginary part accounts for 
absorption. The Laplacian is written as A = (d 2 + d 2 ) for the two-dimensional case which is of interest 
here. The direct solution of this elliptical differential equation is problematic due to numerical instability 
and the fact that the region of interest is large compared to the wavelength (micrometers up to millimeters 
compared to Angstrom). A numerically much more stable alternative is provided by the parabolic wave 
equation, which can be used in controlled and excellent approximation for forward propagation problems 
in the x-ray regime [20121122] , To this end, a function u(x,z) is defined by ^(x,z) —: u(x, z) exp(— ikz). 
Neglecting terms of order 0(d 2 u) which is justified if the z-axis is almost parallel to the wavevector, the 
Helmholtz equation is fullfilled for 4", under the condition that u solves the parabolic wave equation 
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[-2ikd z + di + k J (n 2 - 1)]« = (2) 

Mathematically, this equation is of the same form as the stationary Schrodinger equation for a massive 
particle (without spin) in a potential. After renaming of variables, the potential in the Schrodinger equation 
is equivalent to the the refraction index n profile in the optical case considered here. 

In order to simplify the expressions and to find a suitable length scale for numerical calculations, wc 
use the natural units for the distance along the propagation direction, and W — for the lateral 

dimension, respectively. Note that W is the critical width at which a waveguide becomes mono-modal, i.e. 
the width at which the waveguide supports only a single mode [18 23 . At the same time, a waveguide of 
width d — W leads to the highest possible wave confinement. The intensity distribution of the wave broadens 
both for smaller and larger d. While the second case is trivial, it is important to note that the nearfield 
intensity distribution is broadened by evanescent waves in the cladding, if d is reduced below W. Using these 
natural units 

Z := Skz, X -=^> (3) 



© can be rewritten as 



in material 

d z u = ~dlu+{ 6 (4) 



7T 



0, in vacuum. 



The ratio j remains the only material dependent parameter in this equation. 

The equation is solved using a Crank-Nicolson finite-difference scheme, which was already proved to give 
reliable results for x-ray waveguides and Fresnel zone plates [22124) . For the initial values, a plane wave 
propagating along z was assumed. The lateral boundary conditions are given by a damped plane wave 
propagating in the material far away from the region of interest. 

For field propagation in vacuum, the Fresnel-Kirchhoff integral for two-dimensional beam propagation [25) 
was used [26] . 

i f e~ lkr 
*(x,z) = — / dx'^(x',0) — — cos a , (5) 
V A J V r 

where a is the angle between the z-axis and the vector (x — x 1 ,z). The corresponding distance to the origin 
is r 2 = z 2 + (x — x') 2 . 

3. Tapered x-ray waveguides: shape optimization 

We have considered two generic types of tapering geometry: (i) the linear taper, and (ii) a tapering 
function parameterized by Bezier curves. For comparison with experimental values, we have chosen ma- 
terial parameters and photon energy E according to recent experiments and present fabrication methods. 
Vaccuum channels in Si have been simulated at a photon energy E — 12keV. Taking air/vaccuum as the 
guiding material simplifies simulations because only the index of refraction of the cladding differs from unity. 
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Optically, the structures are still very close to the experiments with polymer channels in Si [7 27J, and are 
completely equivalent to the second generation of x-ray waveguide channels fabricated by e-beam lithogra- 
phy, ion etching and subsequent wafer bonding of a cap wafer. For both tapering geometries, the waveguide 
end width was held constant at d = 1, W = 20 nm which is the critical width W = fc y^ for single-mode 
waveguides [18 23J. At this value, the resolution limit of waveguides is reached, corresponding to a full width 
at half maximum FWHM m j n = 0.64 W for the intensity profile in the channel of width W [18] . 

The goal of the optimization was to quantify the field propagation and confinement, for given open- 
ing width D and end width d or for given d only. In general, one can expect that optimization of x-ray 
nanostructures might be highly non-trivial, as for many optimization problems occurring in physics [28 29 . 
Nevertheless, since we restricted ourself here to the two above mentioned tapering geometries, we could use 
a quite simple optimization algorithm, see below. The two tapering geometries were parameterized as follows 
(see Fig. [H): 

- In the linear case, a straight interface seperates the channel from the cladding at an opening angle The 
length of the waveguide is denoted by L. 

- In the case of a Bezier curve interface, the shape is controlled by two vectors Pi = (xl,yl) and P2 = 
(x2, ?/2), which determine the derivative of the curve at the endpoints and thereby fully define a curve of 
third polynomial degree [30j . 

We have used these functions also because they can be handled by pattern generators of modern e-beam 
lithography units. 

Firstly, let us consider the linear taper. Waveguides with linear taper have been simulated for different 
length L and angle 1?. Fig. [2] shows the resulting peak intensities I normalized to the incident intensity 
Iq, i.e. the gain I/Iq, at the end of the waveguide (exit plane) as a function of L and Hence, this plot 
displays the "gain landscape" of the optimization problem. Even for this two-parameter space, a very rugged 
gain landscape is visible: With increasing opening width D ~ 2Ltan(?9), I first increases, in proportion to 
the geometric acceptance. The initial increase of I{d, L = const) crosses over to a surprisingly broad range 
of angles where the radiation transport and collimation raimains high, but is modulated by an oscillation 
occuring with a periodicity at which small. Finally, the intensity decreases when d becomes too large, leading 
to a leaking of the intensity out of the confining interfaces. This leakage occurs for multireflections of high 
order n. The condition to prevent the leakage is simply given by nd < ?? c , where $ c is the angle of total 
reflection, as expected from ray optics. Nevertheless, at the end of the waveguide a dicrete set of rod-like 
intensity spikes are observed in the simulation, travering or 'leaking' from the interface, see again Fig. [1J 
These rods or spikes originate from interference effects which would correspond to transmitted beams after 
multiple reflections in the ray-optical description. When $ is varied, these rods move and give rise to an 
oscillatory behavior in the exit plane. This effect is characterized by the oscillation between either sharp 
and intense profile at the exit, i. e. high I and small with (full width at half intensity FWHM) or a smeared 
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Fig. 1. Shape functions and field distribution of taperd waveguides: the top row shows the schematics, the bottom row the 
calculated intensity distribution for the optimal parameters of each kind. The waveguides consists of vacuum channels em- 
bedded in Silicon at a photon energy of E = 12keV. (bottom left) Linear taper with an opening angle •& = 0.037 $ c , trans- 
porting radiation collected within a width of D = 0.77/^m over a propagation distance of L = 8 mm down to the critical 
exit cross section of d = 19.88 nm, corresponding to the critical width W of Si. (bottom right) Waveguide tapered with a 
Bcziercurve interface. The lengths d, D and I as above. The two control vectors are P\ = {x\ = 3.45mm,j/i = 766.2 nm) and 
P2 = (x2 = 7.25 mm, yi = 176.7 nm), respectively 

profile with low / and large FWHM. In summary, L = const) possesses an optimum angle z9, while 
/(^ = const, L) increases with the length L, if $ is chosen accordingly, at least over the simulated range in 
L. We chose L = 8 mm as an experimentally achievable upper limit, see Fig. [2] For this length, the optimal 
angle is § ~ 0.037^ c leading to an opening width of D = 0.77/j,m and I/Io — 51.4. 

Now let us consider tapering of the Bezier type. Obviously, more parameters are needed to describe 
the space of Bezier curves which includes the linear as well as the parabolic case. Correspondingly, an 
optimization of the parameters is needed, since the optimum can no longer be found 'by eye' from running a 
set of simulations. As a first step to a complete optimization, we optimized curves by simulating waveguides 
for different y\ and j/2 at fixed x\ — and X2 = \L. As in the linear case the exit was set to the critical 
width d = W. D and L were kept constant. 





0,04 0,08 0,12 0,16 

Fig. 2. (top) Maximum intensity 7(i?, L) at the exit normalized to the incident intensity Iq, as a function of opening angle 
i? and waveguide length L for linear taper, (bottom) Maximum intensity /(■$, L = 8mm) (left axis) and beam width FWHM 
(right axis) at the waveguide exit, as a function of i3. The oscillations are phase shifted with respect to each other. 
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Fig. 3. Maximum intensity I (normalized by Iq) of a waveguide with Bezier-type tapering parameterized by Pi = (l/3L,j/i) 
und P2 = (2/3L,?/2)- The function I(yi,y2) exhibits a well pronounced global maximum. 



The results, see Fig. [3] showed again a rugged gain landscape with a couple of local minima. Nevertheless, 
the maximum intensity seems to occur close to the linear case, and a few local minima are present in 
this particular region. Therefore, a simple approach was feasible: using the amoebae algorithm |31j . the 
optimum parameters Pi, P2 were searched, starting from the solution of the linear case Pi — {^L , ^- — ^A) 
and P2 — (§£, -j — §A) with A = |(D — d). Note that using other starting values, the algorithm also 
found different local maxima, proving the complexity of the underlying gain landscape and the usefullness 
of our pre-scanning of the parameter space. As expected, the maximum intensity is higher for the more 
general Bezier case than for the linear taper, namely 82.2 versus 51.4 for the peak gain I/Io, and F = 72.24 
versus 63.27 for the integrated flux, see Tab. [TJ However, this difference is surprisingly small. The result 
are tabulated in Tab. [TJ showing maximum intensity F, integrated flux J = J dxl, width w (FWHM) of 
intensity in the exit plane, as well as the angular acceptance on the incidence side, defined as the width 
(FWHM) of 1(a), where a is the incidence angle a. Note that if not otherwise stated, the angle a between 
incident beam and waguidc axis was kept constant at a = 0. In addition to the two tapering geometries, the 
planar waveguide is included as a reference. As expected, the angular acceptance u> is much smaller for the 
tapered waveguides due to conservation of photon density in phase space (brilliance). 
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Table 1 

Simulation results: maximum intensity F, flux J, width (FWHM) of intensity at the exit w, and FWHM for the angular 
dependency of the intensity (angular acceptance) ui arc tabulated for (a) the plane, (b) linear tapered and (c) a waveguide with 
interfaces parameterized as Bezier curves (values obtained after parameter optimization). 

(a) planar (b) linear (c) Bezier 



intensity ± 11, 1 • lCT 3 51.4 82,2 

flux ^ 10,4-KT 3 63,27 72,24 

beam width w 15,2 17,8 12,0 

angular acceptance f- 0,929 0,031 0,025 



4. Summary and Conclusion 



We have studied beam propagation and concentation in tapered x-ray waveguides by simulations of finite- 
difference equations and subsequent parameter screening and optimization. The local slope of the waveguide 
interfaces and the associated tapering geometry is limited by the angle of total reflection, which restricts 
the tapering geometry to very elongated shapes. Over a device length on the order of 10mm, an incoming 
beam with a width of 0.5 — 1/im can be concentrated to the minimum width (FWHM) w — 0.65T4 7 ~ 13 nm 
for silicon. The associated expected flux enhancement in the range of 40 — 80 for one-dimensional focusing is 
observed in the simulation, confirming that most of the radiation is indeed transported to the exit aperture. 
For two-dimensional focusing, the gain would be squared, and would thus range up to about 5 x 10 3 , if 
10mm long structures are available. Note that both lithographic wafer processing and hollow glass fibers 
with well controlled shapes are conceivable even on length scales of 100mm. However, let us stay withing a 
more conservative estimate based on a 10mm long waveguide. As a result, the flux density gain, the size of 
the entrance aperture, and the angular acceptance w are fully compatible with an optical scheme, where a 
synchrotron source is imaged by a lense of moderate numerical aperture (mirror system or CRL) onto the 
entrance aperture of a tapered waveguide. Let us consider realistic values for 3rd generation synchrotron 
sources: a source on the order of 100/zm cross section is demagnified by CRL or KB optics of focal length 
/ ~ 1 — 5m positioned 100m behind the source. In this case, the focal spot and the convergence angle 
of the prefocussing optics can be matched to the entrance aperture and the angular acceptance w of the 
tapered waveguide, respectively. In this case the combined gain of the pre-focussing optics and the tapered 
waveguide would lead to an intense and collimated virtual source, which is unlikely to be achieved by a 
single optical device. Most importantly, the simulation results show that the performance of the tapered 
waveguide does not depend on the shape function in a very sensitive way. Indeed, the beam can be funneled 
to the exit aperture surprisingly well already with a linear taper. In sharp contrast to single bounce mirror 
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optics, the tapered waveguide almost forces the convergence of the beam, and is thus much less sensitive 
to figure errors. Nevertheless, we have also observed a rather rugged gain landspace for this optimization 
problem. Since the optimization was only for a few-parameter space, a simple algorithm could be used here. 
We expect that for the optimization of more complex x-ray nanostructures, like layered materials, more 
sophisticated algorithms have to be applied. 
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